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Abstract 

By using the intertwiner and face- vertex correspondence relation, 
we obtain the Bethe ansatz equation of eight vertex model with open 
boundary conditions in the framework of algebraic Bethe ansatz method. 
The open boundary condition under consideration is the general so- 
lution of the reflection equation for eight vertex model with only one 
resctriction on the free parameters of the right side reflecting bound- 
ary matrix. The reflecting boundary matrices used in this paper thus 
may have off-diagonal elements. Our construction can also be used 
for the Bethe ansatz of SOS model with reflection boundaries. 
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1 Introduction 



One of the most important goal of exactly solvable lattice models is to find 
the eigenvalues and eigenvectors of the transfer matrix of a system, then to 
obtain the thermodynamic limit of this system. 

Bethe examined the completely isotropic case of the XXX model and 
found the eigenvalues and eigenvectors of its Hamiltonian[l]. After Bethe's 
work, Yang and Yang analyzed the anisotropic XXZ model by means of 
Bethe ansatz[2]. Then, Baxter in his remarkable papers gave a solution for 
the completely anisotropic XY Z model [3]. He discovered a relation between 
the quantum XY Z model and eight vertex model which is one of the two 
dimensional exactly solvable lattice model. Faddeev and Takhtajan simph- 
fied Baxter's formulae and proposed the quantum inverse scattering method 
or algebraic Bethe ansatz method to solve the six vertex and eight vertex 
models, whose spin chain equivalent are XXZ spin model and XY Z spin 
model, respectively[4] . After Yang-Baxter- Faddeev- Takhtajan's work, a lot 
of exactly solvable models have been solved by algebraic Bethe ansatz [5, 6], 
functional Bethe ansatz[7,8], co-ordinate Bethe ansatz[9], etc. 

Typically, the two-dimensional exactly solvable lattice models are solved 
by imposing periodic boundary conditions in which the Yang-Baxter equation 
provides a sufficient condition for the integrability of the models. 

R\2{ZI - Z2)Rlz{zi - Zi)R2z{z2 - Z^) = R2i{z2 - Zz)Riz{zi - Zi)Ri2{zi - Z2) (1) 

where, the R-matrix is the Boltzmann weight for the vertex models in two 
dimensional statistical mechanics. As usual, Ri2{z), Riz{z) and R2z{z) act in 
C" ® C" ® C" with Ri2{z) = R{z) 0l,R23{z) = 1® R{z), etc. 

The exactly solvable models with non-periodic boundary conditions have 
been early studied in ref. [10-14]. Recently, integrable models with open 
boundary conditions have been attracting a great deal of interests. This 
was initiated by Cherednik[15] and Sklyanin[16], they proposed a systematic 
approach to handle the open boundary condition problems which involves 
the so called reflection equation (RE). 

RuiZl - Z2)Ki{zi)R2liZi + Z2)K2{Z2) 
= K2{Z2)R12{ZI + Z2)Ki{zi)R2l{Zi - Z2) (2) 

The open boundary conditions are determined by the boundary reflecting 
matrix K satisfying the RE. 
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By using a non-trivial generalization of the quantum inverse scattering 
method, Sklyanin obtained the Bethe ansatz equation of six vertex model 
with open boundary conditions by algebraic Bethe ansatz method[16]. The 
transfer matrix with a particular choice of boundary conditions is quantum 
group Uq{sl{2)) invariant [17, 18]. After Sklyanin's pioneering work, a lot of 
exactly solvable lattice models with open boundary conditions have been 
solved. Mezincescu and Nepomechie solved the A[^'' and Ag'^'' vertex models 
by using the fusion procedure[19]. Foerster and Karowski solved splq{2,l) 
invariant Hamiltonian which contains a non-trivial boundary term by us- 
ing the nested algebraic Bethe ansatz [20]. Using the same method, de Vega 
and Gonzalez- Ruiz solved the An vertex model, they also analyzed the ther- 
modynamic limit of this system[21]. Yue, Fan and Hou solved the general 
SUg{n\'m) vertex model[22]. For other progress about open boundary condi- 
tions along this direction, see e.g. ref. [23-31]. 

However, for Baxter's eight vertex model with open-boundary conditions, 
little progress has been made. Jimbo et al obtained the difference equation 
of n point function for semi- infinite XYZ chain [32]. Yu-kui Zhou have 
recently studied the fused eight vertex model and have found the functional 
relations for eight vertex model with open boundary conditions [33], in which 
the reflecting K matrix which satisfy the vertex RE(2) are diagonal. 

In this paper, we will use the algebraic Bethe ansatz method to solve 
the eight vertex model with open boundary conditions. It is known that 
K matrix is a solution of RE. The general solution of RE for six vertex 
model was obtained by de Vega [34]. In ref. [35,36], solutions of RE for eight 
vertex model have been found. The general solution of RE for eight vertex 
model which has three free parameters was found by two groups [37,38]. In 
our approach we use the general solution for left and right boundaries. The 
reflecting boundary K matrices thus may have off diagonal elements. We 
need only to impose one relation on the free parameters of the right K matrix, 
if the free parameters of the left K matrix are arbitrarily given. 

It is known that in Baxter's original work, in stead of Yang-Baxter re- 
lation, the star-triangular relation plays the key role which can be obtained 
from Yang-Baxter relation by using the intertwiners. This was later gener- 
alized to Zn Baxter-Belavin model[39,40] to describe the interaction-round- 
a-face model by Jimbo, Miwa and Okado[41]. This intertwiner method was 
also used to solve the Bethe ansatz problem for similar cases [42,43]. In 
order to get algebraic Bethe ansatz equation of eight vertex model with 



3 



open-boundary conditions, we need to describe the exchange relations of 
the monodromy matrix in the "face language". Thus in our paper, we need 
to convert our boundary conditions of vertex model to that of a face model. 
Our approach is equivalent to a SOS model with open boundaries satisfying 
the face RE. The face RE was first proposed by Behrend, Pearce and Brien 
[44]. In ref. [44], they also find a diangonal solution of face RE for ABF[45] 
model. By using the intertwiners, the face RE is derived directly from vertex 
RE and the general solution of face RE for eight vertex SOS model are found 
by other groups[46,33]. In this paper we actually use a diagonal solution of 
face RE at the left boundary and an upper triangular solution at the right 
boundary for eight vertex SOS model. This open boundary conditions are 
different from the case discussed by Yu-kui Zhou in ref. [33]. Since a diagonal 
matrix is a special case of the upper triangular matrix, our approach can 
be used to a SOS model with boundaries proposed by Behrend et al. We 
can prove that by taking a special case, the Bethe ansatz equation for RSOS 
model, ABE model, can also be obtained. 

The outline of this paper is as follows. In section 2, we first review 
the eight vertex model and reflecting open boundary conditions, the model 
under consideration in this paper will also be constructed. In sect. 3, by 
using the correspondence of face and vertex models, we derive the face RE. 
As mentioned above, the face RE will play a key role in the algebraic Bethe 
ansatz method for eight vertex model instead of vertex RE. The transfer 
matrix with boundary conditions will also be constructed by using the face 
weight. In sect. 4, we will find the local vaccum for eight vertex model with 
a boundary which is the same as the vaccum state given by Baxter[3,7] in 
the Bethe ansatz of eight vertex model with periodic boundary condition. 
The face boundary matrix derived directly from vertex boundary matrix is 
obtained. In sect. 5, the Bethe ansatz problem is solved for eight vertex model 
with open-boundary conditions. Section 6 contains some discussions and the 
further work. 



2 Description of the model 
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2.1 The R matrix 



We first start from tfic R-matrix of tlie eiglit vertex model. Denote a = 
{ai,a2),ai,a2 = 0, 1. Let g and /i be 2 x 2 matrices with elements gu' = 
{-iy6ii>,hii: = 6i+i^i>,i,i' = 0, 1. 

Define 2x2 matrices = I{ai,a2) = h°'^g'^^,Io — I — identity, and 
define I^^ — I ® I ® ■ ■ ■ la® ■ ■ ■ ® I ■,Ia''^^ J-th space. As usual, /^■'^ act in 
V = Vi ® ■ ■ ■ ®Vu where space V is consisted of / two-dimensional spaces. 
We then introduce some notations used in this paper. They are 



e 



a 
b 



h{z) 



^ exp jyT v^(m + a) [(m + a)r + 2(2; + 6)] } , (3) 



m€Z 

e 



1 I Ql 

1 I 0.2 

2 2 



= d 



2 2 
1 



[2,2r),z = 0,l, 



la,(^ + f) 



The R-matrix of eight vertex model takes the form (Fig. la) 



(4) 
(5) 
(6) 

(7) 



(8) 



which satisfy the Yang-Baxter equation (Fig. 2a) 

— ■RjkiZj Zk)Rik{Zi ^kjRijiZi Zj). 



(9) 



It can be proved that the R matrix of eight vertex model satisfy the following 
unitarity and cross-unitarity conditions, 



unitarity : Rij{z)Rji{—z) = p{z) ■ id, (10) 
cross - unitarity : Rl^{z)Rfi{-z - 2w) = p'{z) ■ id, (11) 
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where id is the identity and ti denotes transposition in the i-th space. p{z) 
and p'{z) are scalars satisfying 

p{z)=p{-z), (12) 

p'{z)=p'{-z-2w). (13) 

In eqs.(8-ll), the indices take the value i,j^k — 1, • • • , Z. 



2.2 Reflection equation and reflecting boundary con- 
ditions 

We now deal with an exactly solvable lattice model with reflecting boundary 
conditions, the R-matrix defined above is the Boltzmann weights for this 
lattice model. In order to construct the transfer matrix of this system, we 
must introduce a reflecting boundary matrix K{z) which is a 2 x 2 matrix and 
satisfy the reflection equation proposed by Cherednik[15] and Sklyanin[16] 
(Fig.2b). 



^12(2:1 - ^2)^^1(^1)^21(^1 + ^2)-f^2(^2) 
= K2iz2)Rl2iZl+Z2)Ki{zi)R2liZi-Z2) (14) 

As mentioned in the introduction, two groups have obtained independently 
the general solution K of this reflection equation for eight vertex model. Here 
we take the solution K{z) as [38] 

K{z) = Y.C^^^ (15) 

where Cq, are arbitrary parameters. Correspondingly, we have the dual re- 
flection equation which is necessary in the following of this paper (Fig. 2c): 



Rl2{z2 - Zi)Ki{zi)R2l{-Zi - Z2- 2w)K2{z2) 
= ^2(^2)^^12 (-2:1 - Z2- 2w)Ki{zi)R2l{z2 - Zi). (16) 

We take the solution of this reflection equation as 

K{z)^j:Ca /" (17) 
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where Ca are also arbitrary parameters. Usually, we also call the reflecting 
boundary K and K matrices as right and left boundary matrices, respec- 
tively. 

In order to deal with the systems with open boundary conditions, let us 
define two forms of standard "row-to-row" monodromy matrices >S'i(2;i) and 
Ti{zi) which act in the space y = Vi (8) ^ ® • • • ® V/ by 

S'i(2;i) = Rii{ui + Zi)Ri_i^i{ui_i + zi)- ■ ■R3i{u3 + zi), 

Ti{zi) = Ri3{zi - U3)Ru{zi - W4) • • • Rii{zi - ui), (18) 

where it;_i, • • • M3) = {ui} are arbitrary parameters. S2{z2) and T2{z2) 
can also be similarly defined. 

Considering {-Uj} are the same for Si{zi) and (^2)5 and noticing that 
two R matrices acting on four different spaces commute with each other, we 
find 

i?2i(^i - ^2)<S'i(2;i)5'2(2;2) = ^21(2^1 - Z2)Rn{ui + zi)Ri2{ui ^2) • • • • (19) 
Using Yang-Baxter equation repeatedly, we have 

R2l{Zl - Z2)S^{Z^)S2{Z2) = S2{Z2)S^{Z^)R21{Z^ - Z2) . (20) 

Similarly, we can also obtain 

Tl(Zl)i?12(^l + ^2)^2(^2) = S2{Z2)R12{ZI + 22)^1(^1), (21) 
r2(z2)Ti(2l)i?12(^l - Z2) = Rl2{Zi - Z2)Ti{zi)T2{z2). (22) 

For the periodic boundary condition cases which are studied extensively 
before, the transfer matrix is defined as the trace of the standard "row- 
to-row" monodromy matrix. But for the open boundary conditions cases, 
instead of the standard "row-to-row" monodromy matrix, we should define 
the "double-row" monodromy matrices which take the form: 

h{z,) = T,{zi)Ki{z,)Si{z,), 

k2{z2) = T2{Z2)K2{Z2)S2{Z2). (23) 
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Using the relations listed above, we can prove ki{zi) satisfy the reflection 
equation 

-^12(2^1 - Z2)ki{zi)R2i{zi + 2:2)^2(^2) 

= k2{z2)Rl2{Zl + Z2)ki{Zi)R2l{Zi- Z2). (24) 

As used usually in the framework of the quantum inverse scattering method, 
ki{zi) are 2x2 matrix with elements defined as operators acting in the space 
V = V3 ^ V4 ^ ■ ■ ■ <S) Vi which is the so called quantum space, the spaces 
Vi and V2 are the auxiliary spaces. Eqs.(14) and (24) shows that k{z) is the 
co-module of K{z). 

2.3 The transfer matrix 

Now, let's formulate the transfer matrix with open boundary conditions. 

t{zi) = Try, {Ki{zi)ki{zi)} 

- J2K{z,)kik{z,)ik, (25) 

kl 

with i — 1,2. Since the transfer matrices are defined as the trace over 
the auxiliary spaces Vi,i — 1,2, they should be independent of Vi and V2, 
and are represented as operators acting in the quantum space V3 (8> ■ ■ ■ ® V^. 
With the help of the unitarity, cross-unitarity relations of R matrix, Yang- 
Baxter relation, refiection equation and its dual refiection equation, we can 
prove that the transfer matrices with different spectrum commute with each 
other [47], 

t{zi)t{z2)=t{z2)t{zi). (26) 

This ensures the integrability of this system. 

The aim of this paper is to find the eigenvalues and eigenvectors of the 
transfer matrix which defines the Hamiltonian of the system under consider- 
ation. We will use the algebraic Bethe ansatz method to solve this problem. 
The transfer matrix is defined as a linear function of the elements of the 
"double-row" monodromy matrix. So, it is necessary to find the proper linear 
combinations of the elements of the " double- row" monodromy matrix whose 
commutation relations are suitable for algebraic Bethe ansatz. Besides this. 
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we also need to find an "vaccum" state which is independent of the spectrum 
z. It is well known that this "vaccum" state can be obtained easily for six 
vertex model with periodic- or open-boundary conditions. For eight vertex 
model, it is not a trivial problem. We will study the commutation relations 
and the "vaccum" state problems in the following sections. 

3 Commutation relations 

It is known that for six vertex model and other trigonometric vertex models 
we can obtain the necessary commutation relations directly from the reflec- 
tion equation in which k{z) is the "double-row" monodromy matrix. But 
for the eight vertex model whose R matrix has eight non-zero elements, we 
can not obtain such relations directly from the reflection equation. We have 
to use the vertex-face correspondence to solve this problem. That means 
we should properly combine the elements of k{z) so that we can find simple 
commutation relations which can be dealed with by algebraic Bethe ansatz 
method. 

3.1 Face- vertex correspondence 

We first define a two element column vectors (f)m,ij.{z) , /i — 0, 1, m e Z, whose 



where a = m + 7, 7, G C,k = 0, 1. We call m the face weight, the face 
index, which take values 0, 1. is usually called the three-spin operator. 
It can be proved that we can find row vectors (p, (p satisfying the following 
conditions for generic w, (3, 7. 



A;-th element is [7,41,49] (Fig.lc) 




(27) 




(28) 



where 
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The above relations can also be written in other forms 

1 

<Pm+u,u{z)4)m+i',u{z) — I, 



v=0 



IJ.=0 



(29) 



As usual, I is the 2x2 unit matrix. (Fig.3) 

We define the face Boltzmann weights for the interaction-round-a-face 
model (IRF) as follows[7,41,45,49]: 

h{z + w) 



W(m\z)Z = 



W{m\zy^ 



w{m\zr^: 



h{w) 

h{w{m + 'y) - (-1)''^) 

h{w{m + 7)) 
h{z)h{w{m + 7) — {—lYw 



h{w)h{w{m + -f)) '^^ 

where the face indices fj,, u take the values 0, 1. The other face Boltzmann 
weights are defined as zeroes, so we can see explicitly that for a given face 
weight m, we only have six non-zero face Boltzmann weights at all. Tradi- 
tionally, the face Boltzmann weights for eight vertex SOS model are denoted 
cb h 

Its relation with the notations used in this paper is (Fig. lb): 



as 



d 



m + fi + u 171 + 0' 
m + (i m 



(31) 



The reasons that we use notations (30) are that a lot of zero face Boltzmann 
weights will not appear in our calculation. It is also convenient to compare 
the results of eight vertex model case with the six vertex model case. 

The face Boltzmann weights of IRF model defined above have a relation 
with the R-matrix of eight vertex model which is usually called the face- 
vertex correspondence[3,7,41]. 



^^12(^1 - Z2)(t)m+fL+u,tii.^'L)(t>m+0,v{^2) 
= E ^("^1^1 -^2){Iv/<^I^+^'+f>/,^'(^2)0ji+/i/,;,/(2l) (32) 
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where 0^ denote that it act in z-th space (i = 1,2). 

With the help of the properties of (p, 0, (j) (28,29), we can derive the fol- 
lowing relations from the above face-vertex correspondence relation. 

0m+A,At (^1)^12(^1 - ^2)0^+^,^.(^2) 

= -^2);;[;fi+A/+i>,^/(^i)0l^+/i+i>/y(^2), (33) 

= Y: WHzi - ^2)^'^!;;^^,+,, ,,(^i)S.',.'(^2), (34) 
= j:Wim-il- v'\z, - Z2rX4>^loA^i)€UA^2), (35) 

v'lj.' 

= i: W{m\z^ - ^2)^'<^^i.A.+., ^.(^i)02U',.'(^2). (36) 

All of these relations obtained above have described the correspondence be- 
tween face and vertex models. Usually, we call (p the intertwiner of face- vertex 
correspondence (Fig.4). 

3.2 Commutation relations for elements of the face 
boundary reflecting k matrix 

As mentioned above, in order to obtain a comparatively simple commuta- 
tion relations which can be dealed with by algebraic Bethe ansatz method, 
we should change the "vertex" reflection equations to "face" reflection equa- 
tions, since where the new "R-matrix" - face Boltzmann weights have only 
six non-zero elements instead of eight non-zero elements in "vertex" case. 
For this purpose, by using the three-spin operator (p, we change the vertex 
boundary reflecting matrix to face boundary reflecting matrix, and find the 
commutation relations between the elements of the " face" type monodromy 
matrix, which are useful for the quantum inverse scattering method. 

We first change the matrix K{z) defined in eqs. (16,17) to face boundary 
reflecting matrix i^(m|2;)^, using the unitarity properties of the intertwiner 
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(29), we find 

X [(t)m-ij.+vA^)^m-iX+v,v{^)\} 

V 

where 

k{m\z)l = (j)m,i,{-z)K{z)(l)m-fi+i>A^) (37) 

Thus, the transfer matrix of eight vertex model with open-boundary con- 
ditions can be rewritten as 

t(z) = Tr[k{z)k{z)) 

K ixu ) 

= ^[0m-A+i>,i^(^)^(2;)0m,/.(-^)]^("^k))l 

= Y.kim\zrMm\z);, (38) 

which is true for arbitrary m. Here we also introduce the definition of face 
boundary reflecting k matrix as: 

We call m and m — fi + u the initial and final weight of k{m\z)'^, respectively. 
Thus, we have written out the transfer matrix by using the face form of the 
model. 

Next, we will derive the face reflection equation directly from the ver- 
tex reflection equation by using the intertwiner. Multiply both sides of 
eq.(14) from left by ff+^„,^o(2i)ff+/io+^o,^o(^2), from right by 0^1+^3 ,^3 (-2^1) 
0^+^3+,>3,j/3(— ^2), notice the properties such as 0*^^^ commutes with ^^^^ and 
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ki{zi), use the face- vertex correpondence relations. We can get the face 
reflection equation (see Appendix A and Fig.5,6) [44,46,33], 

k{m + p,2 + i'i\zi)^lk{m + /is + h\z2)ul 

= k{m + fii + i'2\z2Xlk{m + fi3 + i'3\zi)'j^l 

W{m\z, + Z2r,i:iWim\z, - Z2):f^,- (40) 

Here and below summation over repeated indices are assumed. One can find 
that this equation is true for arbitrary ^q,vq, ^^^v^. We should point out 
here that this face reflection equation is different from the one proposed by 
Behrend et al in ref.[44], if the cross- inversion relation of face Boltzmann 
weights [41] is applied, the two equations are equivalent. 

Now, we will let the indices in the above relation take special values so 
that we can find the necessary commutation relations. When Ho = i^o = 
0, A«3 = ^^3 = 1, wc get 

k{m\zi)lk{m - 2\z2)lW{m\zi + Z2folW{m\zi - ^2)00 
= k{m\z2)lk{m-2\zi)lW{m\zi + Z2folW{m\zi- Z2)n. (41) 



1 00 
'00 



From the definition of face Boltznamm weights, we know that W{m\z 
W{'m\z)\\ = -^^j^^y^, so we find 

k{m\zi)lk{m — 2\z2)q — k{m\z2)l,k{m — 2\zi)q. (42) 

This means that the position of Zi and Z2 can be exchanged in this form. 
Let /Xo = i'o = ^^3 = 0, /[/3 = 1, we obtain 

k{m + 2\z2)'^Qk{m\zi)lW{m\zi + Z2)°ooW{m\zi - ^2)^ 
= k{m\zi)lk{m\z2)oW{m\zi + Z2)o\W{m\zi - ^2)00 
—k{m\z2)lk{m\zi)QW{m\zi + Z2)olW{m\zi — Z2) 



Denote 



01 
10 

-k{m\z2)lk{m\zi)lW{m\zi + Z2)llW{m\zi - ^2)^}. (43) 

A{m\z) = k{m\z)Q, 

D(m\z) = k{m\z)l, 

B{m\z) = k{m\z)l, 

C{m\z) = k{m\zfi. (44) 
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So, the matrix ^ ^ D j boundary k matrix in the "face" form. 

Relations (42,43) give the commutation relations of BB and AB, respectively. 
In order to use the algebraic Bethe ansatz method, we must also get the 
commutation relation of DB. 

Let ^0 = H3 = U3 = l,uo = 0, exchange zi and Z2, we get one relation, 
exchange zi and Z2 in equation (43), we find another relation, combine the 
two relations and with the help of equation (43) , we can find the commutaion 
relation of DB: 



k{m + 2\z2)lk{m\zi)lW{m + 2\zi + Z2)o\W{m + 2\zi - Z2)\o 
W{m\z, + Z2)^W{m\zi-Z2yol 
= k{m\z2)lk{m\zifoW{m + 2\zi + Z2)llW{m\zi + Z2)ol 

{W{m + 2\z2- zi)llW{m\z2 - ^i)^^ - W{m + 2\z2 - zi)llW{m\zi - Z2)[l} 
+k{m\z2)lk{m\zi)lW{m + 2\z2 - zi)%W{m\z2 - ^1)0} 

{-W{m + 2\zi + Z2)uW{m\zi + Z2)oo + W{m + 2\zi + Z2)%W{m\zi + ^2)0?} 
+k{m\zi)lk{m\z2)'^oWirn + 2\zi + Z2)%W{m\zi + zs)^!} 

{-W{m + 2\z2 - zi)l\W{m\z2 - 2i)?J + W{m + 2\z2 - zi)%W{m\zi - ^2)^} 

+k{m\zi)lk{m\z2)\W{m + 2\z2 - Zi)\lWim\z2 - ^i)[ji 

{W{m + 2\zi + Z2)nW{m\zi + 22)^0 - W{m + 2\zi + Z2)%W{m\zi + ^2)0?} 

(45) 

In the right hand side of this equation, behind k{m\z2)l, we find not only 
term k{m\zi)\ but also term A;(m|^i)Q, this will cause trouble in the proceed- 
ing of the algebraic Bethe ansatz method, especially in the case where the 
thermodynamic limit of this system is taken. In order to solve this problem, 
we should reformulate A and D as ^ and D so that when we commute D 
with B, only term D exists behind B{z2) which is the notation of k{m\z2)l. 

It is not easy to find D by direct calculation, but the work of Sklyanin[16] 
gives us a hint to formulate the term D. Sklyanin pointed out in his paper 
that the term D in six vertex model case is one of the elements of the inverse 
matrix of the monodromy matrix. Now we will study the inverse of face form 
monodromy matrix for the eight vertex model case. 
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For simplicity, abusing the face weights m etc, we can write the face 
reflection equation as 

WuiZi - Z2)ki{Zi)W2l{Zi + Z2)k2iz2) 
= k2{^2)Wi2{zi + Z2)ki{zi)W2l{zi-Z2). (46) 

It is just the same as the vertex reflection equation. If the inverse of k{z) 
exists, we have 

k2^{z2)Wi2{Zi - Z2)ki{Zi)W2l{Zi + Z2) 

= Wr2{zr + Z2)kr{zr)W2i{zi- Z2)k2\z2). (47) 

That means that k2^{—Z2) and k2{z2) has similar exchange relation with 
ki{zi) in formalism. We have known that the commutation relation of k{m-\- 
212:2)0 with k{m\zi)\ is comparatively simple which is also useful for algebraic 
Bethe ansatz. From the above equation, we assume that the commutation 
relation of k~^{rn + 2|z2)o with k{rn\zi)Q should have the similar properties 
in formalism. We know that A;^^ (2^2)0 is a linear combination of A and D in 
six vertex model. We hope that this is also true for the case of eight vertex 
model. Fortunately, we have obtained the expected result. 
Define 

Q{m\z-wy^^ = k{m + fis + i)2\z-w)';iW{m\2z-wyjl>^lW{m\-w%r, 

Q'imlz-wr.^ = k{m\z-w)>;iW{m\2z-w)'^i:iW{rn\-w)ll^^, (48) 

where, as usual, summation over repeated indices is assumed and = 1,1 = 
0. We can show (Appendix B) 

Q{m\z — w)'^ k{m + u + {^"Iz)"^,' = p{m,i'\z)5^t,f' , 

k{m + p," + il'\z)'^Q'{m\z-wyJ = p'{m, ii\z)5^^>', (49) 

where p{m,u\z), and p'{m,fi\z) are scalars of the "quantum space". We 
can rescale Q and Q' such that p, p' 1. But this is not necessary in the 
following derivation. Multiply (40) by Q{'m + u + jlo\z2 — wYj' from left and 
by Q'{m + /is + v'\z2 — w^,^ from right. Summation over uq, 1/3 gives 

Q(m + z> + po\z2 - wy"k{m + /ig + />2ki)|:jVF(m|zi + Z2yf^l 

xW{m\zi - Z2Y^yj-J{m + As + v\ i^'\z2)5u^u' 

= p{m + z> + /io, v\z2)5y^^k{m + /is + h\ziY^\Q'{m + As + t>'\z2 - w)u3 

W{m\z, + Z2)Z:iW{m\z, - Z2)llZ. (^0) 
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In the derivation of LHS, we use the fact that iy(m|2;i + Z2)'^l'^l 7^ only if 
J^i + IJ^2 = 1^2 + fJ^3- The equation becomes 



Q{m + U + fio\z2 - wYJ'kim + /is + i>'ki);ii 

W{m\z^ + Z2)ZlWim\z,-Z2r,lZ 
= k{m + /is + h\zi)f^Q'{m + /is + t>'\z2 - wY^^ 

W{m\z, + Z2XZ-W{m\z,-Z2)ll% 
p{m + v + fiQ,v\z2) 
p\m + fi^ + v',v'\z2y 

which is similar to (40). Put po = u — u' — 0, fj,^ — 1, and notice m + + = 
m + 2,m + + l = m. We then have 



Q(m + 2\z2- w)lk(m\zifoW(m\zi + Z2)%W(m\zi - Z2)ol 
+Q{m + 2\z2 - wfok{m\zi)lW{m\zi + Z2)olW{m\zi - Z2)[ 

p{m + 2,Q\z2) 



+Q{m + 2\z2 - w)lk{m\zi)\W{m\zi + Z2)olW{m\zi - 2:2)01 



00 
00 
10 



= k{m\z^)lQ'{m\z2 - wY,W{m\z^ + Z2T,^W{m\z^ - Z2Z ' ^, . 

(52) 

Using the same derivation as that in Appendix B, we can calculate 
p(m + 2, 01^2) {-l)h{wa)h{w{a+l)) 



p'(m,0|z2) h{w{a-l))h{w{a + 2)y 



(53) 



where a = m + 7 and 7 is a parameter in defining 0. From (48) and consid- 
ering the explicit form of W , we can write Q, Q' by components of k{z), 

Q{m\z - w)l = k{m - 2\z - w)lW{m\2z - w)llW{m\ - w)°l, (54) 

which is proportional to B{m — 2\z — w), and 

Q{m\z-w)o = -\k{m\z -w)lW{m\2z-w)lo 

-k{m\z - wfoW{m\2z - w)%\ W{m\ - wfol 
= -Q'imlz-wfo 

= -D{m\z-w)W{m\-wyol, (55) 
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where 

D{m\z) = -k{m\z)QW{m\2z + w)il + k{m\z)lW{m\2z + w)lo 
W{m\ - wfol 

is indeed a linear combination of A{m\z) = k{m\z)Q and D(m|2;) = k{m\z)\. 
Subsititute (53-55) into (52) and write k{m.\z2 — w)\ in the equation in terms 
of D{m\z2 — w) and A{m\z2 — w). We then change the parameter Z2 — w io 
Z2 giving the expected equation 

D{m + 2\z2)B{m\zi) 

h{zi - Z2- w)h{zi + Z2 + 2w) 



B(m\zi)D(m\z2] 



+B{m\z2)A{m\zi 



h{zi + Z2 + w)h{zi - Z2) 
^ h(2zi)h(zi + Z2 + w(a + 2))h{2z2 + 2w) 
h{2zi + w)h{zi -\- Z2-\- w)h{w{a + 1)) 



+B(m|..)£)(m|..) "'t; " "t ^ t ""it'"'' ■ (57) 

^ ' ^ ^ ' ^ h{2zi + w)h{zi- Z2)h{w{a + 1)) ^ ' 

where h{z) = ao{z). In the derivation we have used a formular of 9 function 
[7] 

h{u + x)h{u — x)h{v + y)h{v — y) — h{u + y)h{u — y)h{v + x)h{v — x) 
— h{u + v)h{u — v)h{x + y)h{x — y). (58) 

We need also to change D in the exchange relation of A and B by linear 
combination of D and A. Thus (43) is rewritten as 

A{in + 2\z2)B{m\zi) 

h{zi + Z2)h{zi — Z2 + w) 



B{m\zi)A{m\z2) 



—B{m\z2)A{m\zi) 



h{zi + Z2 + w)h{zi - Z2) 
h{2zi)h{zi — Z2 + w{a + l))h{w) 



h{2zi + w)h{zi — Z2)h{w{a + 1)) 



^/ , X -p^/ , N h(—Zi — Z2 + wa)\h(w)]'^ , . 

-B{m\z2)D{m\z^— ^ — , , ^\), '\ -. 59 

^ ' ^ ^ ' ' h{zx + Z2 + w)h{2zx + w)h{w{a^l)) ^ ' 

We have also used (58) in the derivation. With the permutation relations 
of AB, DB and BB (57,59,42), we can obtain the algebraic Bethe ansatz 
equations provided 
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1. The transfer matrix t{z) is a linear combination of A and D, 

2. There is a "vaccum state" of the quantum space which is an eigenstate of 
A and D but not an eigenstate of B. 

We will study these problems in section 4. 

4 Vaccum state and boundary conditions 

4.1 Vaccum state 

The algebraic Bethe ansatz requires to construct a state of the "quantum" 
space , which is an eigenstate of operators A and D with all spectrum pa- 
rameter z. This state is called a vaccum state, which is also an eigenstate 
of D. Before introducing the vaccum state, let us make some preperations. 
We first expresss S and T (see Eq.(18)) by the "face laguage", change their 
auxiliary indices to face indices, and express k{m\z)'j^ by such expressions of 

5 and T. Prom (23,39), the operator k{z) with "face" indices can be written 
as 

k{m\zy^, = 4>m+iu'-i>',n'{z)T{z)K{z)S{z)(f)my{-z). (60) 
From (28,29) we have 

1X1/ 

Combining these two equations gives 

k{m\zy^, = ^K{mo\zy^T{m - z>',mo - i'\z)ij,'ij,S{m,mo\z)t,'i,, (62) 

where we define 

K{mo\zY^ = Aio+f^'oA^)Ki^)^rnoA-z)^ (63) 

T(m - z>',mo - z>|2;)^'^ = (j)m+f,'-u',^^'iz)T{z)(j)„,^^+f,_^^^|J_{z), (64) 

S{m,mo\z)^',, = 4>moA-^)Si^)'Pmyi-z). (65) 

Equation (62) is true for all m and mo- 

Next, assume that we can properly choose the parameters of the right 
reflecting matrix such that for given niQ and all z, K{mo\z)1 = 0. This 
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is possible. We will study this problem in the second part of this section. 
Actually, this requirement constitute the only restriction on the boundary 
matrices in our approach. 

Then, multiply 4>^]i{z)(pinQ+ifl{-z) from left and multiply o+i,o(^)'/'Si(-2^) 
from right to equation (21). Using face- vertex correspondence relations 32- 
36) we can prove the following exchange relation of S and T with face indices, 

W{m + l\2z)\lS{m — 1, mo + l\z)ioT{m, rno\z)io 
+W{m + l\2z)1lS{m + 1, mo + l|2;)oor(m, mo|2;)oo 
= I^(mo + l|2z)°Jr(m + 1, mo + l|2;)ii5'(m, mo|2;)ii 

+W^(mo - l|2^)0°T(m + l,mo - l|^)io^(mmo|^)io, (66) 

which will be useful in our derivation. For later convenience we rewrite S, T 
as 

S{z) = Rio{ui + z)Ri^ifi{ui^i + 2) ■ • • Riq{ui + z), 

T{z) = Roi{z - Ui)Rq2{z - M2) • • • Roi{z - ui). (67) 

They are acting on the space Vq ® Vi ® ■ ■ ■ ® V/, where Vq is the auxiliary 
space. The final preparation is the following observation. In Eqs. (32-36), the 
summation over face indices has at most two terms in eight vertex model. 
In the following cases, due to the non zero condition of W{m\z)'j^^\ there is 
actually only one term in RHS of the equations, 

1. -Eg. (32) when /i — v 

2. Eq. (33) when j-i ^ v 

3. Eq.{3b) when i'- 

Thus we have (Fig. 7) 

-^12(^1 — ^2)</'m+2,o(^l)'/'m+l,o(^2) 
= W{ |^l-^2)|i20SUo(^2)0SUo(^l)> 

'rm— 1,1 Z2)(pm+l,0\Z2) 

= I^(m|zi-Z2)K!i(^i)0So(^2), 



(68) 



(69) 



(70) 
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0m;i(^2)i?i2(^i - ^2)C;o(^i) 

= W^(m|zi-Z2)Si</'SUo(^i)^J^U,i(^2), (71) 

which will be repeatedly used in the proof of vaccum state. 

The above are preparations for introducing the vaccum state. Now, define 
the vaccum state as 

|0 >Z,^ <t>ZM)<t>tAM-^) ■ ■ ■ </'S-(/-2),o(^2)0S_(,_i),o(«l) (72) 

where m — mo — I. This is precisely the same vaccum state introduced 
by Baxter in the original work of Bethe ansatz for eight vertex model with 
periodic boundary conditions [3]. For the vaccum state defined in (72) we 
can show that S'(m, mo|-z)oo and T(m, mo|-z)ii change |0 >™j to |0 >Z~li, 
while S{m,mo\z)ii and T{m,mo\z)oo change |0 to |0 >Z'^li (with some 
coefficients). The operators S{m^mo\z)oi and T{m,mo\z)io change |0 
to zero. Following is the proof. 
We have 



S{m,mo\z)oo 

= <Pml,o{~^)^lo{ui + z)Ri-ifl{ui-i + z) ■ • •i?i,o(Ml + z)(l)^i]o{-zi).{73) 

Since vectors and operators belonging to different spaces may change their 
positions in an equation, we obtain 

5'(m,mo|2;)oo|0 y'Z^ 

= <Pml,oi-z)Rio{ui + 2:)(/)2,,o(«0-Ri-i,o(^^«-i + ^)0mo-i,o(^i-i) 

• • • -R20(«2 + z)(/)^So-{l-2},oMRw{Ui + 2;)0J^+i,o(Mi)</'S!o(-^)- (74) 

By using (69), we move (I)^^]q(—z) towards left step by step. At each step 
we eliminate an R matrix getting a W factor, and change the face weight of 
0*^°^o(~'2) that of ^'^^•'o encountered. Thus we have 

S{m,mo\z)oo\0 >'^^ = • • • -^20(^2 + ^)(/'m+2,o(^^2)0S+i,o(-^)0m!o(Mi) 
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xl[W{ \u, + z)Z (75) 



i=l 



Due to the orthogonal relation (28), this becomes 



i=l 



soo{z)\0 >Z:\ . (76) 



Similarly, one can show 



^(m,mok)oi|0>™ = Cii(-^)Cio(-^)--- 

= 0. (77) 

For the action of Su, from (65,67) and (72), we write 

S{m,mo\z)u\0 >Z, 
= 0mLi(-^)^«o(w/ + z)4)^!l^^f^{ui)Ri-ifl{ui-i + z)</)Si;i\^o('"i-i) 

• • • Rw{ui + ^)0m+l,o(wi)0m!l(-^)- (78) 

Using (71) we move 4>^^\i{—z) towards right step by step. At each step we 
eliminate an R matrix getting a W factor, and change the face weight of 
4>^^\i{—z) and 0*^*^Q encountered. We then have (Fig.8) 



5(m,mo|^)ii|0>- = J]^w{mo-{i-l)\ui + zf^^'^Z-iA-^) 

1=1 

= ilW{mo-{i-l)\u, + z)^,\Q>Z:i, 



i=l 

^ s,M\z)\Q >Ztli. (79) 
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here mo — I = m. Similar derivation can be performed for T, giving 
T{m,m,\zUQ>Z, = \{W{ \z - u,f^\Q >Ztl, 

i=l 

= ^oo(^)|0 >^^+i, 
T{m,mo\zU\Q>Z, = 0, (80) 

T(m, mo|^)ii|0>- = {m + {i - l)\z - Uif^\^ 

i=l 

^ h,{m\^)\0 >Z;1^, (81) 

which completes our proof. 

We now assume that in (62), K{'mo\z)1 = 0, and consider fi = (—1)'^, 
obtaining 

k{m\z)Q = K{mo\z)QT{m - l,mo - l\z)ooS{m,mo\z)oo 
+K{mo\z)lT{m - l,mo + l\z)ooS{m,'mo\z)oi 
+K{mo\z)lT{m - 1, mo + l\z)oiS{m, mo|z)oi, (82) 

k{m\z)l = i^(mo|2;)oT(m+ l,mo - l|2;)ioS'(m,mo|2;)io 
+K{mo\z)lT{m + 1, mo + l\z)ioS{m, 'mo\z)u 
+K{mo\z)lT{m + 1, mo + l\z)nS{m, mo|-2)ii. (83) 
By (66), k{m\z)\ can be written as 

Hm\z)\ = _^("^ok)[i \w{m + l\2z)\\S(m - l,mo + l\zUT(m,mo\zU 
W [mo — i-\^z)QQ >- 

+l^(m + l|22;)?j5(m+ l,mo + l|z)ooT(m,mo|;2)oo} 
+K{mo\z)lT{m+ l,mo + l|z)io5'(m, mo|2;)ii 



K(mo|^)[]W^(mo + l|22; 



01 



xT(m + 1, mo + l|2;)ii5'(m, mo|^;)ii. (84) 

Acting on the vaccum state |0 >Zq and by equations (76,79-81), these two 
operators yield 

k{m\zfQ\0>Z, = K{mo\zf,T{m-l,mo-l\z)ooSoo{z)\0>Z;\+0 + 
= K{mo\z)%,{z)sooiz)\0>Z, 

= r(m|z)°|0>-, (85) 
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and 



k[m\z)^\Q>^^ = 1 [2^)00 «oo(^)*oo(^)|0 +0 

+ {K{m^\z)\- 



K{m^\z)lW{m^ + l\2z)% 



W{m,-l\2zf^ 



xtii(m+ l|2)sii(m|z)|0 



mo 



^ r(m|^)}|0 >;::„. (86) 

From equations (85,86) we see that when K{mQ\z)\ = 0, the vaccum state 
|0 >™j^ is indeed an eigenstate of ^ = k{m\z)Q and D = k{'m\z)\. Thus it 
is also an eigenstate of D. The eigenvalues depend on m, {ui} , I and z, the 
spectrum of A and D. 



4.2 Boundary conditions 

The algebraic Bethe ansatz requires a vaccum state, which needs the right 
boundary to satisfy 

K(mo|2)2 = 0. (87) 

Also, the transfer matrix t{z) must be a linear combination of A and D, 
which impose the left boundary to satisfy 

Kim^Yo = (88) 

and 

K(m^\z)l = 0, (89) 

for mP which we will specify in the next section. We will see that mP is 
constrained with m, mo of vaccum state by 

m + 10 = m + I = mo, 

m+ (-1 + 6)/' = m + 2/' = m°. (90) 

In equation (90), / is the column number of the lattice and I' is a positive inte- 
ger, which will be the number of B operators in constructing the eigenstates 
oit{z) (see (105)). 
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From the general solution of RE (15,17), 



k{z) ^ ^" 



and the definitions 

we can derive the following results (Appendix C), 

^ ' ^ ' h{z + wl3 + w -\)h{wa-w) ^ ' 

- h{-z + wP-'^)h{wa) ' ^^^^ 



^^"^^^^^ - ^ h{-z + wP-\)h{wa) 



(94) 



where a = m + 7. We can easily see from (92-94) that those conditions 
(87-89) are actually independent of z, and depend only on {Cq,} , jC'cvj and 

parameters /?, 7, iw, T,etc.. For any given generic jCaj, we may solve the 
equation 

^C',(-l)"V„(2r7) = 0. (95) 

a 

The LHS of (95) is a doubly quasi-periodic holomorphic function of rj. From 
the quasi-periodicity we see that [7,40,41] it has four zeros in : 77 — >• 
rj -\- l^-q ^ rj -\- T. Assume 771, 772 are two different zeros, we obtain a and /? 
by solving 

-w{a-l) + wj3-'^ = 27/1, 

w{a+l)+wl3-]- = 2772. (96) 
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Then from a = m° + 7 we obtain 7 according to a given m°. Using such l3, 7 
from generic jC^aj, one may construct (p, (p and with which the equations 
(88) and (89) are satisfied. Since (3 and 7 are completely determined, equation 
(87) is a constraint for {C^} at the right boundary. This is the only restriction 
we must impose for our approach. There are 3 free parameters at the right 
boundary for general solution (15), (an overall scalar of K is not important). 
With the constraint (87), there are still 2 free parameters left. If we further 
require 

K{mo\z)l = 0, (97) 
then only one free parameter at right boundary survives. We can show 

K{mQ\z)\ _ h{C - z)h{aw + i + z) 
K{mQ\z)l ~ h{C + z)h{aw + C - z) 

where ^ is the right boundary free parameter introduced in [44] , and similarly 
for the left boundary. In this case the left and right boundary are equivalent 
to that of SOS model (except a symmetric factor) introduced by Behrend, 
Pearce and Brien in [44] where they derive the solutions directly from the 
face reflection equation. This implies that our approach can be used for the 
Bethe ansatz of SOS model with such boundary conditions. 



5 Bethe ansatz 

We see that for a given mP, properly choosing /? and 7, we can ensure 

k{m'^\z)l = K{m°\z)l = (98) 

at the left boundary. We need the parameters {Cq} at the right boundary 
to satisfy 

^a(-l)"V„(u;a + t/;/9-^) =0 (99) 

for an integer mo, where a = mo + 7, to ensure K{mo\z)i — 0. Assume that 
21' = mP — {niQ — I) is a positive even integer, where / is the number of column 
of the lattice. We can proceed the standard algebraic Bethe ansatz [3,4,7,16] 
as following. 
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We have the transfer matrix from Eq.(38), 

t(z) = ^(m°|^)°it(m°|z)° 

+K{m^\z)lk{m^\z)\. (100) 

Due to (44) and (56), t{z) can be rewritten as 

t{z) = tXo{z)A{m^\z) + iii{z)D{mP\z) (101) 

We also have the exchange relations of A, D with B, Eq.(57) and (59). They 
can be compactly written as 

A{m\u)B{m — 2\v) = aoo{m,u,v)B{'m — 2\v)A{m — 2\u) 

+6oo("^) v)B{m — 2\u)A{m — 2\v) 
+6oi(m, v)B{m - 2\u)D{m -2\v), (102) 

D{m\u)B{m — 2\v) — aii{'m,u,v)B{'m — 2\v)D{m — 2\u) 

+bio{m, u, v)B{m — 2\u)A{'m — 2\v) 
+bn{m, u, v)B{m - 2\u)D{m - 2\v). (103) 

Besides, we have (42), the exchange relation of 5's, 

B{m\u)B{m - 2\v) = B{m\v)B{m - 2\u). (104) 

Consider a vector of the quantum space 

$ = S(m° - 2|2i)5(m° - A\z2) ■ ■ ■ B{m' - 2l'\zi,)\0 >Z,, (105) 

where the number of B is V , rmP — 21' — m. Due to (104), it is symmetric 
for Zj's. We will show that for properly chosen zi, - ■ ■ , zi' which satisfy the 
so called Bethe ansatz equations, $ is an eigenvector (eigenstate) of t{z). 
We have 

t{z)^ 

= {fXo{z)A{m!'\z) + //i(^)L>(m°|z)} $ 
= {fXo{z)A{m^\z) + Hi{z)D{m^\z)} 



mo 



xB{m° - 2\zi) ■ ■ ■ B{m\zi>)\0 > 
{S(m° -2|zi) \lio{z)aQQA{mP — 2\z) + fj,i{z)a[iD{'mP — 2\z) 
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+S(m° - 2\z) [^o(^)feoo^(^° - 2|^i) + ^o(^)&oi^("i° - 2^0 

+//i(^)6;o^K - 2|^i) + /Xi(^)6'nL'(m° - 2\z^)\ } 

xSK-4|z2)---|0>;^„ 
= {B(m° - 2|zi) [(//a')o^(m° - 2\z) + (iia')iD{m'^ - 2\z)] 

+B{m^ - 2\z) [{fib')oA{m^ - 2\zi) + {^b')iD{m'^ - 2\zi)] } 

XSK-4|Z2)---|0>™ 
= {S(m° - 2\zi)B{m° - 4\z2) [(/xa'a")o^(m° - 4\z) 

+ {lia'a")iD{m^ -4\z) 

+B{m^ - 2\z)B{m^ - A\z2) [{^Jih' a'[\A{m^ - A\z{) 
+{lib'a'l)D{m^ - A\zi) 

+B{m^ - 2\z)B{m^ - A\zi) [{^b'b")oA{m^ - 4\z2) 
+ ifib'b")iD{m'> - A\Z2)]} 



x5(m° -6|^3)---|0 > 



m 
mo' 



(106) 



where a', a", a'/, 6', 6" are 2x2 matrices. The matrices a',a",a'( are diag- 
onal. The elements of these matrices are determined by mP, z and Zi via 
Eqs.(102,103). The notations of these matrices are 

a{mP,z,Zi) = a', 
a{mP — 2, z, Z2) = a", 



a{m!^ — 2,zi,Z2) = 

b{mP,z,Zi) = b', 
b{m^-2,z,Z2) = b" 



1) 



(107) 



for short. The notations (// • • ■)i represents the z-th component of the product 
of the row vector fj, with matrix In the following derivation, the 

notations are similar. Repeatedly using Eqs.(102,103) to move A and D to 
the right of all S's, we obtain 

t(z)^ = B(m'^ -2\zi)---B{m'^ -2l'\zi>) 

X Uiia'a" ■ ■ ■ a(''))oA(m° - 2l'\z) + {i^a'a" ■ ■ ■ a^^\D{rnP - 2l'\z) 
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+S(m° - 2\z)B{m^ - 4\z2) ■ ■ ■ B{m^ - 2l'\zv) 



X 



1° 



{lib'a'i---a^P)oA{m'' -2l'\zi) 

+{fj,b'a'( ■ ■ ■ a{\b{wP - 2l'\zi)\ |0 >l 
+B{m'^ - 2|z)B(m° - 4|zi)S(m° - 6|^3) • • • 

A(m° - 2^1-^2) + • • • - 2l'\z2)] |0 + 



The vector |0 is an eigenvector of A{m\z) and D{m\z). The eigenvalues 
are given by (85,86). Thus we have 



A{m\z)\0>Z, = Ao(^)|0>-, 
D{m\z)\0>^^ = Ai(^)|0>- . 



(109) 



Therefore, noticing m — 21' = m, we can write t{z)^ as 
t{z) = Uiia'a"---aP)oXo(z) + (i^a'---a^%X,{z) 



m 
mo 



xS(m° - 2|zi) • • • B{m\zi>)\0 > 
+ [ilJ,b'a'i---af\Xoizi) 

+ ifib'a';■■■aP)^X^iz^) 
xB{m^ - 2|z)5(m° - A\z2) ■ ■ ■ B{m\zv)\Q >\ 
+ [(• • •)oAo(^2) + (■ ■ ■)iAi(^2)] 5(m° - 2|^)S(mO - A\z^) 



mo 



xB{m'-6\zs)---\0>Zo + 



(110) 



Because of (102,103) and (104) we see that t{z)^ must be a hnear combina- 
tion of 



B(m''-2\z,)---B{m\z,)\0>^^ = 



B{m^ - 2\z)B{m° - ^z^) ■■■\0> 



mo 
m 
mo 



*1 



i?(m°-2|^)i?(m°-4|zi)---i^(---ki-i)i?(---ki+i)---|0>;;:o = * 



5(m° - 2\z)B{m^ - 4\zi) ■ ■ ■ S(m|^,/_i)|0 > 



m 
mo 
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This is because we can always change the order of S's such that ZiS inside 
5's are arranged according to the order of i. So we have 

t{z)<^ = C^^o + Cl^i + Cl^2 + --- + Cl,^i' (112) 

The problem now is that although the forms of Cq and Cl are simple and 
clear, for i > 2 are represented by a complicated summation. However, 
using the fact that $ is a symmetric function of z'^s, we can greatly simplify 
the calculation. Let us exchange Zi and Zi in $. This does not change $. 
Then we can use the above standard procedure to have 

t{z)^ = C^* + Ci^i + • • • + Ci^i + + • • • (113) 

where Cq and Q can be obtained by exchanging Zi and Zi in Cq and C}. 

Assume ^ii ■ ■ ■ ) are linearly independent vectors. Then each coeffi- 
cient for the linear decomposition of t{z)^ by {^^j} is unique. Thus we have 
Cl = Cl Put all = for i ^ in (112) we have 

t{z)^ = Cl^ = t(z)$. (114) 

It is, $ is an eigenstate of the transfer matrix t{z) with eigenvalue 

T{z)=i^{z)a'a"---a^''^X{z) (115) 

The spectrum parameters Zi are determined by the I' conditions — 0,i — 
1, - ■ ■ ,1'. The first condition is 

d = ii{z)b'a'l ■ ■ ■ afhizi) = 0. (116) 

Other /' — 1 conditions can be obtained by exchanging zi and Zi in (116). 
These are the Bethe ansatz equations. Using the explicit form of (101-103) 
(i.e. Eqs.(56), (100), (57) and (59)), we can prove that these equations are 
actually independent of the spectrum parameter z. This implies that $ is an 
eigenstate of all transfer matrices with arbitrary spectrum. 

To end this section, we present here some results. The left boundary 
matrix is diagonal, we can show that its diagonal elements can be explicitly 
written as: 

K{m°\zfo = h{{a^ - l)tv)h{^- z- w)h{{a^ + l)w + ^+ z)F{z), 
K{m°\z)l = h{{a^ + l)w)h{z + i+w)h{{a° -l)w + i- z)F{z), 
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where ^ is the left boundary free parameter. We notice that this solution is 
identified with the solution given in ref.[44]. While F{z) is a function of z 
depending on the scale of the left boundary matrix, which is not essential. 
So we get 

^o(^) = h{2z + 2w)h{^- z)h{z + ^ + a'^^v)F{z), 
Hi{z) =h{w) h{z + i + w)h{i- z+ {a^ -l)w)F{z). 

The eigenvalue of the transfer matrix of eight vertex model with open bound- 
ary conditions is 

Sj h{zi + z)h{zi - z + w) 
t=i H^i + z + w)h[zi - z) 
, I w I \'r( ^'(^i ^ ^ ^ w)h{zi + z + 2w) 

+lll{z)Xiiz) — — r . 

Here {zi} should satisfy the Bethe ansatz equations: 

h{zi + Zj)h{zi - Zj + w) 



[117) 



^l(^i) j=l jj^i H^i ~ ^3 ~ '^)K^i + + 2^y) 

h{w)h{i - Zi + w{a^ - l))h{i + Zi + w) 



for i — 1,- 

Xo{z) = K{mo\z)^o]l 



h{2zi)h{i + Zi + wa^)h{i - Zi) 
, I'. From (85,86) and the definitions of A, D, we have 

' ^h{z + Ui + w)h{z — Ui + w) 



(118) 



X^{z) = 



i=l 



[h{w)Y 

o h{2z + w{ao + l))h{w) 
h(w(ao + l))h(2z + w) _ 

h{z + Ui)h{z — Ui) 



K{mo\z)\-K{mo\z)l 

I 



h{2z + w)h{w{aQ + 1)) 
h{w)h{wao) 



where a° = mP + 7 = + 21', gq = uiq + 7 = + /, a = m + 
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6 Discussions 



We can obtain our trigonometric limit as the following. The intertwiner (or 
three spin operator) defined in (27) is 



1 _ k 

2 2 
1 



{z + {-l)^w{m + 7) + w(3, 2r) 



(119) 



Define 7' = 7+ 2^,/5' = /3+2^ and note /i = 0, 1. Equation (27) reads 



2u; ' ^ 1-^ ' 2w 
1 k 



2 2 
1 



where 



z + (-l)^t/;(m + 7') + wp' + (/u - l)r, 2r) 
{z + {-lYw{m + i) + wP',2T), (120) 



- 27ri ( )( r+z;+ (- 1)^ i« (m+7' ) +«)/5'+ 5 ) 



(121) 



is independent of k. We then rescale (j) to (f)' — i'^cj). At the same time, 
we must perform a gauge transformation which change W ot W' to en- 
sure the face-vertex correspondence. When r ioo , (p'^ ^{z) S^k- W 
goes to a trigonometric R matrix, which is different with the R matrix in 
ref.[16] only by a constant factor. The boundary condition K{m\z)'l — 0, 
K{rn\z)\ — K{m\z)Q — (if we add K{m\z)\ — 0) also approach that of 
ref.[16]. Thus we can show that the trigonometric limit of our model is that 
of ref.[16]. It is a six vertex model with integrablc reflection boundaries. In 
such model, the number of i?'s (/') in the Bethe ansatz state $ is arbitrary. It 
is reasonable that the eigenstate of transfer matrix with maximum absolute 
value of eigenvalue (ETMM) is a Bethe ansatz state. Each such six vertex 
model can be attained as a limit of a sequence of eight vertex model with 
reflection boundaries. In this limit procedure, the Bethe ansatz equations, 
vaccum states and operators A,B,C,D are all approaching that of six vertex 
model. Thus we can reasonably assume in this sequence of eight vertex mod- 
els there is a sequence of Bethe ansatz states which approach the ETMM of 
the six vertex model. Thus it is quite possible that these Bethe ansatz states 
are ETMM of the eight vertex models, especially when they are very close 
to the limit six vertex model since the eigenvalues of the transfer matrix t{z) 
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are discrete [7]. The fact whether a Bethe ansatz state is with the maximum 
absolute value of eigenvalue, should not depend on continuous parameters if 
there is no phase transition in the procedure, also since eigenvalues are dis- 
crete. For the above six vertex model, when / is given (the column number 
is given), the true discrete varible is I'. The above sequence of Bethe ansatz 
states should have same Thus, it is reasonable that in our approach, 
when the left and right boundary condition determine a proper the Bethe 
ansatz state has a eigenvalue of maximum absolute value, which is the most 
important state in thermodynamics. 

We know that Johnson, Krinsky and McCoy calculated the energy of 
excitations of the XYZ model after Baxter obtains the Bethe ansatz of eight 
vertex model with periodic boundary conditions [48]. Now, using the results 
presented in this paper, we may also calculate the energy of the excitations 
of the XYZ model with boundaries. By analyzing the eigenvalues of the 
transfer matrix, one may also get the boundary free energy, thermodynamic 
limit and finite size corrections. Other physical phenomena are also worth 
of studying suh as surface critical exponents and scaling, the central charges 
in conformal field theory etc.. It is well known that eight vertex model is 
equivalent to SOS model, if one impose some restrictions on SOS model, we 
can obtain the restricted SOS model (ABF model). So, if we impose some 
restrictions on the Bethe ansatz of eight vertex model with boundaries, we 
should find the Bethe ansatz for the ABF model with boundary conditions. 
All of these work worth be studied in the future. 
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Appendix A Quadratic relation of the com- 
ponents of the face type k{z) 

The left hand side of equation (24) is Rr2(zi — Z2)ki{zi)R2i{zi + Z2)k2{z2). 
We multiply it from left by 0S+Ao,Mo(^i)'/'2+/io+i>o,^o(^2) and multiply it from 
right by 0j^^^3,^3 (-^1)01^)^^3+^3,^3 (-Z2) obtaining 

LHS = 0;;U,,„(^i)0l^U+^o,.o(^2)i?i2(^i - ^2) • • • . (122) 
Using (34) we can eliminate R12 to have 

LHS = M^(m|^i-Z2);:;-0W^,+,,,^,(^i)^2U,.i(^2)A;i(^i) 

^21(^1 + ^2)^2(^2)</'!^+/i3,^3(-^l)<^lJ+A3+i>3,^3(-^2). (123) 

Move ^(2) over ki{zi) and move over A;2(-22), so that they are at the left 
and right side of R21 respectively. LHS becomes 

LHS = ■■■0l^U,i^i(^2)i?2i(2;i + ^2)</'!i+A3,;,3(-2;l)••• 
= • • • w{m\z, + ^2)::?;:^<^^U+^i,M.(-^i)^SU+A2,..(^2) • • • . 

(124) 

Since the Boltzmann weight W is non zero only if i>2 + fis — 1^1 + 1^2, we have 

LHS^-.. ct>^X.+,,,J-^i)Wim\z^ + Z2):r,l^^Xs+o,,'^.i^2) ■■■■ (125) 

Now ki{zi) and k2{z2) all have their "own" and at left and right side. 
By definition one conclude 

LHS = k(m + fi2 + i>i\zi)';ik{m + fi3 + h\z2)'Ct 

xWim\z, - Z2r,i:iW{m\z, + Z2)':^. (126) 

The derivation of RHS is similar. 

Appendix B Left and Right inverse matrices 

of k{z) 

When n = 2, we have 

Ri2{z) = R2i{z) = Ri2{zy'''. (127) 
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When z — —w^ we have 

P_(12)i?i2(-«;) = i?i2(-H^-(12) = Ri2{-w), (128) 

where P_(12) is the anti-symmetric operator of space Vi ® V2 satisfying 
P_(12)^ = P_(12). On the other hand the R matrix have the property 

ili^i^Rr3{^\ifr'\il^'r' = ^.(^)- (129) 

Thus 

jfP_(12)i?32(^2)P3l(^l)P-(12)[/f]-^ 

= p_(i2)/(:^)P32(^2)[/^)]-^/fi?3i(^i)[/(:^)]-^p-(i2) 

= P_(12)[/(^)]-lp32(.2)/f [/«]-^P3l(^l)/«P-(12) 

= P_(12)[/(^)]-M/«]-^P32(;^2)i?3i(^i)/(^)il^^P-(12). (130) 
It is not difficult to show 

il^)7^)p_(12) = (-1)-+-P_(12) = P-(12)[/^rM/^^]-'- (131) 

Thus 

/(3)P_(12)P32(^2)i?3l(^l)i^-(12)[/F]-' 
= P_(12)P32(-Z2)P3l(-Zl)P-(12) = C/. (132) 

The operator U acting on V^i (8) V2 ® V3 is invariant under the similar transfor- 
mation by which implies that U is equivalent to a unit operator on the V3. 
From (127) we see that this is also true for U' = P_(12)P23(^2)Pi3(^i)P-(12). 
We then consider the RE (14) for the case zi = z — w, Z2 = z, and have 

G — Ri2{—w)ki{z — w)R2i{2z — w)k2{z) 

= Ri2(-w)Ti(z - w)Ki(z - w)Si(z - w)R2i{2z - w)T2{z)K2{z)S2{z) 

Due to (21), 

G = Ri2{-w)Ti{z ~ w)Ki{z - w)T2{z)R2i{2z - w)Si{z - w)K2{z)S2{z) 
= Ri2{-w)Ti{z - w)T2{z)Ki{z - w)R2i{2z - w)K2{z)Si{z - w)S2{z). 

(133) 
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Rewrite Ri2{—w) as P^{12)Ri2{—w), and move Rui—w) towards the right. 
Each time when it goes over a pair of RiiR2i in 21 T2 by YBE (9), we rewrite 
Ri2{—w) as Pl{12)Ri2{—w) and leave Pi. Then we obtain 

G = [P_{l2)R2^Rr,P_{l2)R2^Ri^P-{l2)---P_{l2)R2iRnP-{l2)] 
Ri2i-w)K^{z - w)R2i{2z - w)K2{z)S,{z - w)S2{z) 
= [M]x---. (134) 

Move Ri2{—w) to the left side of Si by RE. One has 

G = [M]K2{z)Ri2{2z - w)Ki{z - w)R2i{-w)Si{z - w)S2{z). (135) 

Similarly, we move R2i{—w) step by step to the right side of S'iS'2 obtaining 

G = [M]---[N]R2i{-w) 
[N] = [P_{12)Ri2RnP-il2)P_{12)Ri_i,2Ri-i,iP-{i'^)--- 

P_(12)i?32i?3ii'-(12)]. (136) 

From the property of U and U', we see that G is proportional to an identity 
operator in the "quantum" space V' — V3 ® ■ ■ ■ ® VJ. Using the similar 
derivation as in Appendix A, we multiply 0(^)0^^) from left and multiply 
0{i)(^(2) ^Yom right of G, and conclude that when zi = z — w, Z2 — z, both 
LHS and RHS of (40) are proportional to identity operator in the quantum 
space V\ Properly choosing indices and noticing M^(m| —w)^^ = ■ ■ ■ S^pSf^'p', 
we have (49). 

Prom the above derivation, we see also that G is anti-symmetric to the 
classical (auxiliary) indices of space Vi<S>V2 and is independent of m. Thus 

p'{m,i^\z) = 4>in+^fii^ - w)4>l+p+f,Jz)Gtf 

xC+i,i(-^ + ^)'^'I+i+o,o(-^) 

= {4>l{z - w)4>'!..{z) - 4>''..{z - w)4>l{z)} 

{<l>l(-z + w)<f>'.(-z) - cfP.i-z + w)<i>l{-z)] (137) 

Similarly, p(m, v\z) can also be expressed as G^ multiplied by a factor which 
depends only on 0, 0. Therefore the ratio is independent of G^. It 

is completely determined by 0,0. Direct calculation gives (53). 
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Appendix C Derivation of the boundary 
condition 



Prom definition we have 



2 2 
1 



{z + {-lYwa + w(5,2T) 



e^^\z + w{{-lYa + 13)) 

e^^\z + x). 



(138) 



where a = m + 7. When cti, 0:2 are integers, by the expression of 6 function, 
we have 



(139) 



giving 



K^'g^^ct^mA^) = (-l)"^e2'^^(^)(^^+^+^+5+"2)0^^^(^ + ceiT + a2). (140) 

On the other hand, from the property of zeros of doublely quasi-periodic 
holomorphic funcion, one can show[49] 



Det 



^(0)(^l) ^(0)(Z2) 

e^^\z^) e('\z2) 



r-i I, / ^1 + ^2 — 1 N , / — 2:2 - 

C X h\ )h\ 



(141) 



where C is independent of z\ and zi- If we write 



(142) 



then %n-vp.,iX'^ elements of its inverse matrix. Thus for any column vector 
^f, the quantity B = J2k ^m+i li^)^'' can be written as 



1 

DetA 



Det 



^m+6,o( 



'z) 

^1 



(143) 
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Combining (140) and (143) gives 



Det 



{z + wa + w(3) 6'(°) {-z + wa + w(3 + aiT + cts) 
{z + wa + w(3) (9(^) {-z + wa + wp + aiT + as) 

(144) 



Met A 

Two determines can be obtained from (141), thus 

aa{wa + wjS — i)(TQ,(— 2;)(— 1)"2 



B' 



[h{z + wl3 + w — \)h{wa — w)] 



(145) 



Substituting the definition oi K{m\z)1 and the expression of ^(z) (15,63), we 
have 



K{m\zf, = (-1) 



E„C„(-l)°V«(«;a + t«/3-i; 



h{z + wf] + w — ^)h{uia — w) 

Other elements of K{m\z) can be similarly obtained. They are. 
K{m\z)Q = a — > —a in RHS of the above equation. 

1 



(146) 



K{m\z)l = 



h{z + wP + w — \)h{wa — w) 

^a2^Q(^/^ + w - \)aa{-z + wa-w) 



(147) 



K{m\z)\ — a 
K{m\z)l - 



-a in RHS of equation (147). 
1 



h{-z + wP - |)/i(wa) ^ 
K{m\z)l = a — > — a in RHS of equation (148). 



J2 C„(-l)°V„(-wa + w + wP--) 



(148) 



/i(-z + W/S - |)/?.(wa) 
K{m\z)\ = a — > —a in RHS of equation (149). 



aa{wP - l)aaiz + wa) 
aa{z + w) 



(149) 
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